Inelastic effects on the transport properties of alkanethiols 
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Using first-principles approaches we investigate local heating and the inelastic contribution to 
the current for various alkanethiols sandwiched between metal electrodes. In the absence of good 
heat dissipation into the bulk electrodes, we find that the local temperature of the alkanethiols 
is relatively insensitive to their length. This is due to the rates of heating and cooling processes 
scaling similarly with length. On the other hand, when considering heat dissipation into the bulk 
electrodes, the local temperature of alkanethiols decreases as their length increases. We also find 
that the inelastic scattering profile displays an odd-even effect with length which compares well with 
experimental results. This effect is due to the alternating direction of the CH3 group motion with 
respect to current flow with increasing C atoms in the chain, and is very sensitive to the structure of 
the carbon-sulfur-gold bond. Inelastic scattering profiles can therefore help illuminate the bonding 
configuration of molecules to metallic surfaces. 



There is an ever-increasing interest in charge trans- 
port in organic molecules due to their potential applica- 
tion in electronic devices. H H 01 II H Recently, 
several laboratories have reported consistent data in the 
resistance of alkyl chains. H H 0, 0, Some the- 
oretical calculations support these results, [13J at least 
in the shape of the current-voltage (I-V) characteris- 
tics. 0, 0| This seems to suggest that a reproducible 
contact can be created between the alkanethiols and the 



electrodes. H E! El El llil ll& However, sev- 
eral current-induced mechanical effects such as forces on 
ions |2£} and local heating 0, [51 13 can generate sub- 
stantial structural instabilities which can lead to atomic 
geometries quite different than those assumed theoreti- 
cally. 

In this paper, we focus on one current-induced me- 
chanical effect, namely, inelastic scattering in alkancthiol 
molecular junctions. Our intent is to explore (i) the de- 
pendence of local heating on the length of the alkyl chains 
and gain insight into their stability under current flow, 
and (ii) determine the inelastic contribution to the cur- 
rent. We find that the local temperature of alkyl chains 
is smaller the longer the chain provided that there is 
good thermal dissipation into the bulk electrodes. This 
is due to the insulating character of the alkanethiols and 
is in contrast to results for metallic quantum-point con- 
tacts. |H [2l| We also find that the inelastic scattering 
profile displays an odd-even effect which is very sensitive 
to the structure of the carbon-sulfur-gold bond. This 
odd-even effect is in agreement with high resolution elec- 
tron energy loss spectroscopy (HREELS) experiments on 
the same systems. j24| This demonstrates that inelastic 
spectroscopy can be quite effective in determining the 
atomic-scale geometry of molecular junctions. |25| 

We investigate these inelastic effects in the framework 
of static density-functional theory (DFT) in the local 
density approximation. |2fij The calculations proceed as 
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FIG. 1: Local temperature as a function of bias for various 
lengths of alkanethiols. (a) No heat dissipation into the bulk 
electrodes is taken into account, (b) Heat dissipation into the 
bulk electrodes is taken into account. 



follows: first, the stationary scattering wavefunctions of 
ethanethiolate sandwiched between two bulk gold elec- 
trodes (represented with ideal metals, jellium model, r s s» 
3) are calculated by solving the Lippmann-Schwinger 
equation self-consistently. [26j Due to their insulating 
character, the current in these systems decreases expo- 
nentially with length d as I = I exp(-fld). [HE! El 
By exploiting the periodicity in (CH 2 )2 groups of the 
alkyl chains, we can then calculate the wavefunctions of 
the different molecules by a simple scaling argument. The 
vibrational mode energies and the transformation matrix 
which contains the character of the modes (longitudinal 
versus transverse with respect to current flow) [22] are 
evaluated using total-energy calculations. [3] Through- 
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FIG. 2: Differential conductance (left axis) and absolute value 
of d 2 I/dV 2 (right axis) as a function of bias for various alka- 
nethiols with different numbers of carbon atoms. Due to 
the large relative magnitude of the inelastic features of low- 
energy modes with respect to the high-energy modes, the 
curve d 2 I/dV 2 has been multiplied by a factor of ten for bi- 
ases larger than 50 mV. 

out the paper, the angle that the S-C bond makes with 
the surface normal has been fixed at 43 degrees, which 
was obtained by relaxing the structures at zero bias. This 
angle is in reasonable agreement with the one found in 
previous theoretical work. |29j Below we will show the 
sensitivity of the inelastic current to the variation of this 
angle. The electron-phonon coupling constant for each 
atom and mode of the chain is finally calculated as re- 
ported in Ref. |22j. With these quantities we can then 
evaluate both the local temperature of the junction and 
the inelastic I-V curve. 

Local heating - Let us first investigate the effect of 
heating due to electron-phonon interaction. This interac- 
tion allows for the exchange of energy between electrons 
and the lattice via absorption and emission of vibrational 
modes. Details of the theory of local heating in nanoscale 
structures can be found in Refs. l22l.l3cl Here we just men- 



tion that there are two major inelastic channels that lead 
to a given local temperature in a nanojunction. One is 
due to inelastic processes that occur in the atomic region 
of the junction, i.e., involving the atoms of the alkanethi- 
ols and few atomic layers of the bulk electrodes. Elec- 
trons incident from the right or left electrode can absorb 
(cooling) or emit (heating) energy because of electron- 
vibration scattering with the ions of the junction. The 
other channel is due to dissipation of energy into the bulk 
electrodes via elastic phonon scattering. In quasi-ballistic 
systems, where the inelastic electron mean free path is 
large compared to the dimensions of the junction, elastic 
phonon scattering is the most efficient way of dissipating 
heat into the bulk electrodes. [i2,|3(J Let us first focus on 
the inelastic scattering contribution assuming heat dissi- 
pation into the electrodes is negligible. This can be the 
result of, e.g., weak coupling of vibrational modes local- 
ized in the junction with the continuum of modes of the 
bulk electrodes. pHLE^ 

Denoting by the power absorbed (emitted) 

by electrons incident from the left (right) via a vibra- 
tional mode ^, the total thermal power generated in the 
junction can be written as the sum over all vibrational 
modes of the above four scattering processes |25|: 

This power can be expressed in terms of the electron- 
vibration coupling in the presence of current. |2q When 
the heating processes {W^ 2 and W^' 2 ) balance the cool- 
ing processes (W^' 1 and W^' 1 ), i.e., P = 0, a steady- 
state local temperature is established in the junction. 
This temperature is plotted in Fig. ^a) for various alka- 
nethiols of different length assuming zero background 
temperature. It is clear from the figure that the local 
temperature without heat dissipation depends weakly on 
the length of the molecule. This can be understood quite 
easily since the rate of energy transfer W"' k scales sim- 
ilarly with length for all modes and processes (heating 
and cooling). In other words, even though the power 
per mode decreases exponentially with length by a factor 
exp(— (3d) due to the corresponding exponential decrease 
of the current, the total power (Eq. ^) will be zero at 
approximately the same temperature, as the above expo- 
nential term factors out. 

We now allow for the energy stored locally in the junc- 
tion to dissipate away into the electrodes via the cou- 
pling between the normal modes of the molecule and the 
phonons of the electrodes. As in our previous work, [2^ 
we estimate this thermal conductance assuming the junc- 
tion forms a weak mechanical link with a given stiffness 
K. |3lJ The thermal current into the electrodes via elastic 
phonon scattering is thus 

hh = -j— / deeN L (e) N R (e) [n L {e) - n R {e)\ , (2) 
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FIG. 3: Differential conductance (left axis) and absolute value 
of d 2 I/dV 2 (right axis) as a function of bias for undecanethi- 
olate. The curve d 2 //d^ 2 has been multiplied by a factor of 
ten for biases larger than 50 mV. 



where n^rm is the Bose-Einstein distribution function 
and Nl(r) (e) is the spectral density of phonon states at 
the left (right) electrode surface. Since the stiffness K is 
proportional to e? _1 , [HI the thermal current Ith is pro- 
portional to dr 2 . Therefore, while the local power in the 
junction (Eq.^) that needs to be dissipated decreases ex- 
ponentially with length, the thermal power into the elec- 
trodes that allows for this dissipation (Eq. |2J) decreases 
only algebraically. This implies a lower temperature in 
the junction as a function of length [see Fig. QJb)]. We 
note that this trend for alkanethiols is opposite to that 
found in atomic wires [20L l2l| and is due to their insulat- 
ing character (conversely the current in atomic wires is 
relatively independent of length [13 )■ This result com- 
bined with the fact that current-induced forces decrease 
with increasing length [gjj, [U HH suggests that longer 
alkanethiols are more stable against current flow. This 
is consistent with recent experimental results on similar 
systems. || 

Inelastic current - The coupling between electrons and 
molecular vibrations will cause small discontinuities in 
the conductance at biases corresponding to the energy of 
normal modes. 0, H^, 0] The position and magnitude 
of these discontinuities varies significantly from mode to 
mode and can be measured in experiments such as in- 
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FIG. 4: Absolute value of A 2 I/dV 2 as a function of bias for 
two different surface-S-C angles of undecanethiolate. Top: the 
structure with an angle of 43 degrees. Bottom: the structure 
with an angle of zero degrees. The change in angle transforms 
the mode at 176 mV (indicated by an asterisk) from quasi- 
transverse (top-left schematic) to quasi-longitudinal (bottom- 
left schematic) so that it appears as a new feature in the 
inelastic conductance. 



elastic electron tunneling spectroscopy (IETS) 0, Ilif 
and HREELS. 0,0] The theory of inelastic scattering 
as we use it in this paper is reported in Ref. HjJ We 
just recall here that the relative strength of inelastic fea- 
tures is mainly determined by the character of the modes: 
those that have large longitudinal character with respect 
to current flow contribute most to the features in the 
inelastic conductance. 0, [2f| Q 

In Fig. |2 we plot the inelastic contribution to the con- 
ductance and its first derivative as a function of bias for 
various alkyl chains. A broadening of 1 meV correspond- 
ing to elastic phonon scattering has been introduced to 
make the curve d 2 //dV 2 finite. We show the modes that 
contribute the major peaks in Fig. [3] for the case of un- 
decanethiolate. The major vibrational modes in Fig. [21 
are in good agreement with the results of HREELS ex- 
periments on self-assembled monolayers. [l9|,|2^ The dif- 
ferences between the reported inelastic curves and the 
theoretical ones can be partly ascribed to modification of 
vibrational modes due to molecule-molecule interaction, 
an effect which we do not consider in this work. Most of 
the peaks at small bias (V < 50 mV) are associated with 
stretching, bending, and twisting motions of the alkane 
backbone. The largest peak at low bi as ( also observed 
in IETS and HREELS experiments [301) corres Ponds 
to a rigid motion of the whole molecule with respect to 
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the electrodes. We also observe a weak odd-even effect 
in the spectrum of d 2 I/dV 2 at wavenumbers 850 cm -1 , 
1107 cm" 1 , and 1420 cm" 1 . This effect is due to the al- 
ternating direction of the CH 3 group motion with respect 
to current flow with increasing C atoms in the chain. It 
is most pronounced for the 1420 cm -1 mode, which is in 
agreement with HREELS experiments. [24| 

To conclude we discuss how inelastic spectroscopy can 
probe the atomic-scale geometry of the alkanethiol-gold 
bond. In the above, we have considered a geometry where 
the S-C bond forms an angle of 43 degrees with the sur- 
face normal (see top- left schematic of Fig. 0J. However, 
recent calculations have found different energetically sta- 
ble structures depending on the details of the computa- 
tional methods employed. [11111111113,111 The actual 
geometry will be determined by experimental conditions, 
including the way the contacts are formed. In Fig. 0] 
we show that a change of the above angle from 43 de- 
grees to, e.g., zero degrees can change substantially the 
character of the modes in the structure. In this case, the 
change in angle transforms the quasi-transverse character 
of the mode at 1420 cm -1 to a quasi-longitudinal one (see 
Fig. 0J. This mode is the one responsible for the most 
pronounced odd-even effect we discussed above. There- 
fore, with a surface-S-C bond angle of zero degrees, this 
effect would be reversed as a function of the number of 
C atoms in the chain. Inelastic spectroscopy, combined 
with this type of detailed calculations, can therefore 
help illuminate the bonding configuration of molecules 
in metal-molecule-metal structures. 
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